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ABSTRACT 



This thesis is concerned xvith proving that the eigen- 
values of a specific unsymmetric matrix are real and positive, 
without actually computing them. The method of finite 
differences is applied to the vibration analysis of a 
cantilever beam and leads to an unsymmetric stiffness matrix 
in the eigenvalue problem formulation. The technique 
employed in the proof is based on a perturbation theory 
given by Wilkinson for real symmetric matrices. Application 
of the theory is made to the cantilever beam eigenvalue 
problem. The results verify that the eigenvalues of this 
and other \msymmetric matrices can be proven real and 
positive without their actual values being calculated. 

A multiple-variable formulation of the cantilever 
beam vibration analysis is also examined to illustrate that 
it implicitly contains all of the properties of a symmetric 
matrix. 
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I . INTRODUCTION 



The vibration analysis of many structures can be 
accomplished by the application of finite differences to 
the partial differential equations governing the free 
imdamped vibration of the structure. Assuming the motion 
of the structure to be harm.onic, with natural frequency CO 3 
leads to the matrix equation for the eigenvalue problem 

2 

K0= CD M(^ (1-1) 

vjhere K and M are square matrices and are referred to herein 
as the stiffness matrix and the mass matrix respectively. 

H o f* c Tia-l-nvcl V-Thon K ”nd M P..V& 

symmetric and K or M is positive definite^ the eigenvalues 
2 

CD j are real; when both K and M are positive definite^ 

CD^ is always positive |^Ref. ij . VJhen K and M are developed 

using energy techniques, they are always symmetric; however, 

when they are developed for the same problem by applying 

finite differences to the governing partial differential 

equations of motion and the boundary conditions, they 

frequently emerge numerically as unsymmetric, calling into 

2 

question vjhether numerical results for CD will in fact be 
real. 

In this thesis an uns 3 nmnetric stiffness matrix is 
developed by finite differences for the cantilever beam. 

A perturbation theory for real symmetric matrices given 
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by Wilkinson j^Ref. 2j is used to establish that all of 
the eigenvalues of the unsymmetric matrix formulation are 
real and positive; hence the results obtained for the 
natural frequencies are real, representing undamped harmonic 
motion as assumedo 
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II. DERIVATION OF THE MATRIX EQUATIONS 



The structure selected for analysis is the cantilever 



beam of length L shown in Figure 1. 




FIGURE 1 



The equation governing the free undamped vibration of the 
beam can be given in the form 



in V7hich t is time, and 




( 2 - 1 ) 



$= Y/L 



X = X/L 
ja = p/(Er^) 

where Y is the deflection of the beam, X is the axial 
coordinate, p is the mass density, E is the elastic modulus, 
and r is the radius of gyration of the cross-section. 
Assuming harmonic motion. 



= cp> sin ojt 



( 2 - 2 ) 
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Substituting Equation 2-2 into Equation 2-1 leads to 

-7-% ( 2 - 3 ) 

d X 

^ 2 
in which X=fj.(jJ . 

The boundary conditions on the beam are 

® at X = 0 (2-4a) 

and 

= d^cb _ 0 

dx2 dx3 at X = 1 (2-4b) 

Let the beam be divided into ten equal intervals of 
length hL as sho\m in Figure 1, The nodes j’=-l, j=ll^ and 
J=12 are fictitious nodes located off of the beam. Let 
the derivatives in Equation 2-3 and Equations 2-4 be 
approx'i.mauec bv 



/d±\ ; ( .0 + 0 ) / (2h) 

y<ax j. ^ ^ J+i 



(2-5a) 



1^^ 1 ( _ 20, (2-5b) 



\dx^ A 



(^] . 

W3A 



J , 



(-0 + 0 . -0 +0 ^ ) ( 2 - 5 c ) 

^ %--2 ^ j -1 0 ]- f-l J -+2 



= (0 -40 -h6 0 -40 -H 0 )/h^ (2-5d) 

\dx^/. J"2 j-1 j j+1 j-f2 



in vfhich h=0.1. The truncation error of each of the finite 

2 

difference approximations in Equations 2-5 is 0(h ). 
Applying Equation 2-3 at j=l,2,3^ 10^ Equations 
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2-4a at j=Oj and Equations 2-4t at j=10; and eliminating 
the variables 
the matrix equation 

’ 7 -4 1 0 0 

-4 6 -4 1 0 .... 0 

1-46-410. . .0 

0 1 -4 6 -4 1 0 .. 0 

0 01-46-410.0 

0 .01-4 6-4 100 

0 ., 01 - 46 - 410 . 

0 ... 0 1 -4 6 -4 1 

0 .... 0 1 -4 5 -2 

[ 0 0 2-42 

where I is the identity matrix 
matrix in Equation 2-6 is not symmetric. However, if the 
last rovr of Equation 2-6 is divided by tvra, the resulting 
formulation is symmetric. This results in the mass matrix 
M no longer being equal to I, since the diagonal element 
in the last rovj is equal to one-half; hovrever, M is still 
symmetric and positive definite, and all eigenvalues are 
real. The symmetric stiffness matrix vjill be referred to 
hereafter as Note that all of the finite difference 

equations used to construct are centered about j, 

^It can be shovm that this symmetric formulation 
is identical to a matrix equation developed using energy 
techniques with one-half of a mass interval at j=10. 



9^2 eventually leads to 




Note that the stiffness 



13 



An alternate matrix formulation of the governing 
difference equations can be constructed in a similar manner 
by using the O(h^) approximation j^Ref. 3 




= 0+^20 -10?zS+ 3?6 l/(2h3) 

J ^ j-2 0"1 j J+1 



(2-7) 



instead of Equation 2.5c. 

As a consequence of the fact that Equation 2-7 is not 
centered about j, the fictitious node J=12 is not required, 

and Equation 2-3 is applied at J=l,2,3j 9» Tbe 

boundary conditions are satisfied by applying Equations 2-4a 
at J=0 and Equations 2 -4b at J=10. Eliminating 

and 0 -^^ leads to the imsymraetric matrix formulation 



7 

-4 

1 

0 

0 

0 

0 

0 

0 



-4 

6 

-4 

1 

0 



1 

-4 

6 

-4 

1 

0 



0 

1 

-4 

6 

-4 

1 

0 



0 

1 

-4 

6 

-4 

1 

0 



. . . 0 

0 

0 . . 0 

10.0 
-4100 
6 -4 1 0 

- 46-41 

1 -3^ -l\ 

0 i -1 i 




to. 



Xh'*! 



h 



0 



( 2 - 8 ) 



The stiffness matrix in Equation 2-8 vjill be referred to 

as K . 

2 
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III. ANALYSIS AND RESULTS 



A perturbation theory given by Wilkinson j^Ref. 2j for 
real symmetric matrices will be used to prove that the 
eigenvalues of the real unsymmetric formulation given by 
Equation 2-8 are real. The brief description of the theory 
which follows is expanded in the Appendix. 

Let 

K = A + B 
2 



where is the real unsymmetric square matrix of order n, 

A is an n X n real symmetric matrix, and B is an n x n real 
unsymmetric matrix. Every eigenvalue of lies v;ithin at 
least one of the n circular discs whose centers are the 



c: ven va lue.-. 02 



X ‘Ji CXr'cu j_o.x vxj SCS 



|jB[| 2 ^ the 2-norm of B, where 



H 



B|| 2 = (maximum eigenvalue of B B) 



-H 



and B is the Hermitian conjugate matrix of B. The matrix 
H 

B B is then Hermitian, the eigenvalues of which are all 
real. If any m of the n circular discs form a connected 
domain isolated from the others, there are precisely m 
eigenvalues in this connected domain. In the case of real 
unsymmetric matrices , any complex eigenvalues must occur 
as conjugate pairs. If the i^^ circular disc is isolated 
and contains only one eigenvalue, this eigenvalue must 
therefore be real. 

The matrix can be decomposed in many ways to obtain 
the matrices A and B. One technique is to use the symmetric 
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matrix as A and obtain B by subtraction. This method 
was tried and resulted in such a large ||b|| ^ that only 
two eigenvalues of could be proven real. Consequently 
two other approaches vrere considered in an attempt to 
develop a smaller norm size that v;ould prove all the eigen- 
values to be real and positive. These two approaches are 
referred to as Method I and Method II. 

In Method I, the matrix is decomposed such that B 
contains the smallest number of elements possible, i.e.. 



. 0.-4 6 

...2 -4 

. . . 0 1 



-4 

6 

-3 



2 0 

-4 1 




0 



-1 



2 



-4 6-4 1 0 

1-4 6-4 1 

0 1-4 4|- -1 

. . 0 1-1 I 



K 

2 



+ 



. . . 0 


0 


0 


0 


. . . 0 


1 

U 


0 


-^4 

1 


. . . 0 


1 

2 


0 


Oj 



B 



I 



Equation 3-1 shows only that portion of that requires 
alteration. The eigenvalues of and the 2-norm of B^ 
were determined using the IBM SSP subroutine EIGENP and 
are tabulated in Table I and plotted in Figure 2. Note 
in Figure 2 that five of the circular discs are isolated 
and in the positive domain. As a result, applying the 
perturbation theory proves that at least five of the nine 
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eigenvalues are real and positive. Also note that the 
five real eigenvalues are the five largest eigenvalues . In 
order to prove that the lower eigenvalues in the overlapping 
circular discs are real, the inverse formulation 



1 I ^ = k“^ 0 

X 2^ 



( 3 - 2 ) 



is examined. If 1 /^ can be proven real and positive, then 
is real and positive. 



TABLE I 

EIGENVALUES AND N 0 R 14 S - METHOD I 



i 


Eigenvalue \ : 


At 


K2 


At' 


K2-I 


T 

: 1 


-0. 07 Q 

1 - - ; 1 


0 . non s 

1 ■ ' 1 


0.067 

1 ■ 1 


0.066 

1 : 


2 


0.222 


0.049 


0.083 


0.079 


3 


0.361 


0.369 


0.119 


0.107 


4 


1.19 


1.31 


0.209 


0.168 


5 


3.10 


3.16 


0.472 


0.316 


6 


5.97 


5.96 


1.57 


0.764 


7 


9 . 4 o 


9.36 


10.63 


2.71 


8 


12.72 


12 . 68 


91. 4 i 


20.42 


9 


15.12 


15.11 


720.2 


800.1 



Value of 
Value of 



B. 



B- 



=o. 8 i 4 



=180.9 



17 



FIGURE 2. 

METHOD I RESULTS 'N THE COMPLEX -PLANE 
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EIGii; LVALUES OF 



METHOD I RE<;iJLTS IN THE COMPLEX-PLANE 
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AlsOj inversion directly relates the lowest eigenvalue of 

-1 



K 



and so on. Since the 



Kg to the largest eigenvalue of 

lower eigenvalues of K2 have not yet been proven real, the 

higher eigenvalues of Yi ^ are of particular interest. 

The decomposition of K“^ for Method I is 

2 



1 

2 


1 




2 


pi 

‘-■ 2 . 


3 


3 # 


4 


7 


1 


3 


5 


7 


9 


11 


13 


15 


26| 


li 


5 


9 l 


14 


m 


23 


27 # 


32 


57 


2 


7 


l 4 


22 


30 


38 


46 


54 


97 


2 | 


9 


18| 


30 


42 # 


55 


0 ( 2 


80 


i 45 


3 


11 


23 


38 


55 


73 


91 


109 


I 99 i 




13 


27? 


46 


67# 


91 


Ii 5 i 


l 4 o 


259 


4 


15 


32 


54 


80 


109 


l 4 o 


172 


322 


4i 


17 


35 # 


kjcL 


92 | 


12Y 


16^: t 


l:u4 


3or 



K 

2 



1 

2 


1 


• 






• 


• 


4 


4 # 


1 


• 


• 








• 


15 


17 


• 




• 








• 


32 


36# 


• 


• 


• 








• 


54 


62 


• 


• 


• 








• 


80 


92# 


• 


• 


• 








• 


109 


127 


• 


# 


• 








• 


l 40 


164 - 


4 


15 


32 


54 


80 


109 


l 4 o 


172 


204 


hi 


17 


36 # 


62 


92# 


127 


164 # 


204 


U) 

00 
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( 3 - 3 ) 



0 


• 


0 


• 








0 




0 


• 




• 








0 


9 k 


0 


• 




• 








0 


20 ^ 


0 


• 




• 








0 


35 


0 


• 




• 








0 


52 i-| 


0 


• 




• 








0 


73 k 


0 


• 




0 








0 


94 | 


0 


• 




• 








0 


118 


0 


• 




• 








• 


0 



The eigenvalues of A i and the 2-norm of B , are also 

I I 

shown in Table I and in Figure 3 . Examination of Figure 3 

-1 

reveals that the largest eigenvalue of is real and 

positive; hence, the lowest eigenvalue of Kg is real and 
positive. Thus, six of the nine eigenvalues have been 
proven real and positive. To verify that these eigenvalues 
are in fact real and positive, the actual eigenvalues of 
and were computed using EIGENP and are tabulated 
in Table I and plotted in Figure 2 and Figure 3 . They do 
in fact lie within the circular discs. 

The decomposition of K 2 and of Method I led to a 

relatively large norm and the proof that only six of the 
nine eigenvalues vjere real and positive. Consider nov7 
Method II in which is decomposed such that the moduli 
of the elements of B are as small as possible (making 
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B skev;-symmetric ), 



i.e 



. 6-410 

. -4 6 -4 1 

. 2 - 3 % 4| -i3/4 

. 0 i -1 i 



(3-4) 



K 



2 



6-410 

• ® % 

. 1 H -i3/8 

• ° \ i 



+ 



0 0 0 0 
0 0 -1/g i 

0 Vg 0 

0 -i % 0 



B 

II 



Only that portion of requiring alteration is shown in 

Equation 3-4. Using the same procedures as in Method I 

leads to the results tabulated in Table II and plotted in 

Figure 4. The norm of B is much smaller than the norm 

II 

of B^_, and examination of Figure 4 reveals that seven of the 
nine eigenvalues of are real and positive. Decomposing 
in a manner similar to Equation 3-4 leads to 4^^' 
and I . The results are also shovm in Table II and in 
Figure 5* Using this formulation, the lowest eigenvalue 
of Kg is proven real and positive. Consequently, use of 
Method II proves eight of the nine eigenvalues to be real 
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and positive. The last eigenvalue must also be real since 
it can have no conjugate. It must also be positive since it 
is larger than the lowest eigenvalue which has been proven 
positive. Therefore all the eigenvalues of K 2 are real and 
positive . 

A third method considered is to treat in a manner 
similar to the conversion of Equation 2-6 to the symmetric 
formulation K^. In this method^ is altered by multiplying 
row 9 by two. Although this leads to another unsymmetric 



matrix, a much smaller |B ^ be developed. However, 

from past experience it is concluded that this approach 
would not prove all eigenvalues real and positive without 
resorting to the inversion of the matrix. Since Method II 
-hbp piftenvalues of K real and positive without 
alteration of the matrix, and since the eigenvalues of this 
third matrix formulation are not jjjio j this approach is not 
pursued further. 

As a check on the accuracy of the calculated eigenvalues 
of a comparison of the exact natural frequencies Ref. 4 
and the numerical natural frequencies is made in Table III 
for the five lo’west frequencies. 
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TABLE II 



EIGENVALUES AND NORMS - METHOD II 





Eigenvalue \\ 


i 


^II 


K2 


A ' 
^‘II 


Kg-l 


1 


-0.144 


0.0013 


-5.58 


0.066 


2 


o.o4o 


0.049 


0.070 


0.079 


3 


0.363 


0.369 


0.084 


0.107 


4 


1.31 


1.31 


0.119 


0.168 


5 


3ol7 


3.16 


0.209 


0,316 


6 


5.97 


5.96 


0.494 


0.764 


• f ! 


1 9 ♦ 37 


i r 1 


“1 ^0 

! -»- • f ^ 1 


0 Tl 

! - • « - i 


8 


12.68 


12.68 


17. 10 


20.42 


9 


15.11 


15.11 


810.6 

1 


800.1 



Value of 
Value of 



B 



II 



B 



II 



= 0.148 



= 90.3 
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METHOD n RESULTS I'J THE COMPLEX-PLANE 
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I E1GE>JVALUES OF 
X EIGENVALUES OF 
llBlI = 0.148 



FIGURE 5 

METHOD U RESULTS \U\ THE COMPLEX-PLANE 
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TABLE III 



NATURAL FREQUENCIES OF CANTILEVER BEAM 



Natural Frequency rad/sec 


Exact |Ref. 4 j 


Numerical 


0.560 


0.562 


3.510 


3.522 


9.820 


9.671 


19.20 


18.210 


31.80 

i 


28.300 

1 _ . .. -I 
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IV. MULTIPLE VARIABLE FORMULATION 



An alternate unsyrametric formulation is obtained when 
the two variables 0 and "7?^ ^ where is the bending moment 
in the beam_, are used instead of only 0 . The governing 
equations are 






> 2 
c) X 



and 






(4-1) 



(4-2) 



The boundary conditions for the two-variable formulation 
are 



c6 



d0 



= 0 



t-c V/ y V — ' ^ 



and 



d^ 

dx 



=0 



at x=l 



( 'i - -s i 



(4-4) 



Consider the beam shovm in Figure 1. Applying Equation 

4-1 and Equation 4-2 at J=l,2,3,.. 10, Equation 4-3 

at j=0. Equation 4-4 at J=10, and eliminating 0^^ 0^,* 
^ll^^-l^^o^ and^Q_-]_ leads to the matrix equation 
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-2 0 10 . 

-2 -h^ 1 0 0. 

0 1 0-201 

1 0 -2 -h^ 1 0 

0 0 0 1 0 -2 

2 

0 0 1 0 -2 -h 

. 0 0 

. 0 1 



1 0 0 0 0 , 

0 0 0 0 0 . 

0 0 1 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 1 0 0 

0 0 0 0 0 0 0 

I 

0 0 

• • 

. . 0 0 



(4-5) 
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01 \ 
















• • 












02 


• • 












^2 


1 


0 . 


• 








03 


0 


0 . 


♦ 










0 


-2 


0 


1 


0 


0 




• / 

08 


-2 


-h2 


1 


0 


0 


0 




00 


0 


1 


0 


-2 


0 


1 




09 1 


1 


0 


-2 


-h® 


1 


0 




^-^9 1 


• • 


0 


0 


2 


0 


0 




010 


o 


0 


0 


0 


0 


1 




\^ioj 


• 












h\ 


• • 














• • 












02 


• • 


• 










W 2 ' 


0 0 


• 


• 








/ 


0 0 


• 


• • 






( 


\ 

\ • / 


1 


0 


0 


0 


0 


0 




08 


0 


0 


0 


0 


0 


0 






0 


0 


1 


0 


0 


0 


1 


1 

1 


0 


0 


0 


0 


0 


0 


1 

1 


^9 


• • • 


0 


0 


0 


1 


0 




010 



o] \ ^ 10/ 



• 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 



0 0 0 



0 



Equation 4-5 is unsymmetric . However, if the variables 

y?^2> ^10 eliminated. Equation 4-5 reduces 

to Equation 2-6. As previously discussed. Equation 2-6 is 
simply a symmetric matrix in disguise, so that the eigenvalues 
are real and positive. 

The two-variable formulation has been presented here to 
point out that many unsymmetric matrices may in fact be 
reducible to symmetric matrices j hence, they have all of 
the properties of symmetric matrices. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The basic concept illustrated in this thesis is that 
for a given unsymmetric matrix the eigenvalues may be 
examined for proof that they are real and positive vxithout 
actually calculating them. Use of the perturbation theory 
given by Nilkinson gives a relatively straightf or^^rard method 
by Vfhich this may be accomplished. For the simple example 
selected calculation of the eigenvalues is easily carried out, 
thus confirming the results. 

It is recommended that this technique be appD.ied to other 
eigenvalue problem formulations of more involved physical 
systems. If the techniavie proves itself reliable in appli- 
cation to matrices of large order, significant computer time 
can be saved by ascertaining whether a given m.odelling method 
will result in real, positive eigenvalues prior to conducting 
lengthy eigenvalue computations. 
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APPENDIX 



PERTURBATION THEORY GIVEN BY WILKINSON 



With slight changes of notation and text, the following 
material is reproduced for ease of reference from the pages 
cited for Reference 2. 

Consider an eigenvalue \ (possibly complex) of the real 
system 

(A + B)x = .Xx (A-Ia) 

in which A is symmetric but B is not. Since A is real 
there is an orthogonal matrix H such that 

h‘\h = diag (X^ ) (A-lb) 

and theX^ are real. Because X is an eigenvalue of (A+B), 
the matrix (A+B-Xl) is singular and hence its determinant 
is zero. From Equation A- la and Equation A-lb 

h"^(A+B-Xi)H = diag (X^- X )+h'^BH (A-2) 

and taking determinants, the matrix on the right of Equation 
A-2 must also be singular. Two cases are distinguished. 



Case 1. X = Xj_ some i 

Case 2. X?^X-5 any i, so that 



diag(X^- X)+h"^BH = diag(Xi-X ) [ I+diag(Xj_-X)”''H'"BH J (A-3) 

and taking determinants again, the matrix in brackets must 
be singular. If (I +Z ), vrhere Z is a general matrix, is 



-1 -1 



singular. 



A 1; for if I Z 2< none of the 
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eigenvalues (I + Z ) can be zero 



Hence 



Ref. 5 



diag (X^-X)"V^3H 



(A-4) 



giving 



i.e. (in both Case 1 and Case 2) 



max lx.” Xl 


H 

1 




bI 




H 


1 1 '^1 




2 


1 


2 





> 



1 (A-5) 



min 


x.-x 


4 


1 — 1 
1 

tx: 




H 




B 










2 




2 





(A-6) 

for at least one for every eigenvalue \ . That is to 
say^ every eigenvalue of (A + B) lies in at least one of the 
circular discs given by 

-1 



X 



H 



H 



B 



(A-7 ) 



and since H is orthogonal (all eigenvalues of which have 



rilC d'Cl i V S ^ "h Vi P <5 o T* P f1 "n p P T. O 



X.” X 


:< 


B 


1 







(A-8) 



Further^ if any s of these discs form a connected 

domain isolated from the others^ then there are precisely 

s eigenvalues in this domain. 

In the case of a real unsymmetric matrix such as (A + B), 

any complex eigenvalues must occur in conjugate pairs. If 

th 






> 2 



2 (J 7 ^ i), then the i circular disc is 
isolated and contains only one eigenvalue. This eigenvalue 
must therefore be real. Hence, if all eigenvalues of A 



are separated by more than 2 
(A + B) are real. 



B 



,, all eigenvalues of 
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